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Abstract 

We calculate the nuclear dependence of direct photon production in hadron- 
nucleus collisions. In terms of a multiple scattering picture, we factorize the 
cross section for direct photon production into calculable short-distance par- 
tonic parts times multiparton correlation functions in nuclei. We present 
the hadron-nucleus cross section as A"' times the hadron-nucleon cross sec- 
tion. Using information on the multiparton correlation functions extracted 
from photon-nucleus experiments, we compute the value of a as a function 
of transverse momentum of the direct photon. We also compare our results 
with recent data from Fermilab experiment E706. 
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I. INTRODUCTION 



As early as in 1970s, it was observed [|T| that inclusive cross sections for single high 
transverse momentum particle produced in hadron-nucleus scattering show an "anomalous" 
nuclear dependence, in which the cross section at fixed transverse momentum grows ap- 
proximately as y4" with A as the atomic number of the nuclear target. The value of a is a 
function of transverse momentum, and can be as large as 4/3. This phenomenon has been 
known as Cronin effect. Typical energy exchange in a high transverse-momentum scattering 
process is so large that any single hard scattering should be very localized within a single 
nucleon. Consequently, a linear A dependence is expected for single scattering processes. 
Therefore, the Cronin effect is often described as due to multiple scattering of partons in 
nuclear matter The A'^^^ behavior signals a dependence on nuclear size, and multiple 

scattering is dominated by double scattering. 

In some of the previous work on this topic, an independent scattering picture was adopted 
0. In this picture, each scattering was treated independently. For example, the cross section 
for double scattering was proportional to a product of two Born cross sections. The double 
scattering cross section in this picture is not infrared safe however. This is because the 
kinematics of single particle inclusive cross section can only provide a constraint on the 
total momentum from the target, which leaves the possibility that one of the Born cross 
section diverges when the momentum transfer of this Born cross section approaches to zero. 
Therefore, theoretical predictions from the independent scattering picture are sensitive to 
the infrared cutoffs which must be introduced in the calculations. 

Recently, Luo, Qiu and Sterman (LQS) have shown that the anomalous nuclear enhance- 
ment can be described naturally in perturbative QCD, in terms of a nonleading power, or 
"higher twist" formalism 0. In this treatment, contribution from double scattering can be 
factorized into short-distance hard parts convoluted with corresponding multi-parton matrix 
elements, or multiparton correlation functions in nuclei. The short-distance partonic parts 
are calculable in perturbative QCD, and all infrared divergences associated with soft rescat- 
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terings in perturbation theory can systematically be absorbed into multiparton correlation 
functions. The multiparton correlation functions are nonperturbative, just like the parton 
distributions in the single scattering processes. These correlation functions in nuclei provide 
information about nuclear matter and its interaction with high energy probes. They can 
reveal information different from what normal parton distributions in nuclei can provide; 
and in principle, they are as fundamental as the parton distributions. In order to test the 
theory, we need to find different processes which depend on the same multiparton correla- 
tion functions. Information on these new correlation functions extracted from one set of 
processes may be applied in other processes. 

Our aim in this paper is to show that the consistent perturbative QCD treatment of 
double scattering developed by LQS can be naturally applied to high transverse momen- 
tum direct photon production in hadron-nucleus scattering. We factorize the cross section 
of direct photon production into calculable short-distance partonic parts times multiparton 
correlation functions, which are the same as those derived in Ref. ||^. We calculated the 
short- distance partonic hard parts. We evaluate the nuclear dependence by using the infor- 
mation on the multiparton correlation functions, extracted from experiments on momentum 
imbalance of two-jet photoproduction on nuclear targets 0. Our numerical results are con- 
sistent with recent measurement of nuclear dependence in direct photon production from 
Fermilab E706 experiment 0. 

A double scattering with high momentum transfer must have at least one hard scattering 
to produce the high transverse momentum observables. In addition, there may be a soft 
scattering either before or after the hard scattering (referred to below as a soft-hard process), 
or another hard scattering (called as a double hard process). We shall show that only the 
soft-hard processes contribute to the nuclear dependence of direct photon production to the 
order we consider. The fact that the photon does not interact strongly once produced at 
the hard collisions eliminates final-state multiple scattering in direct photon cross sections. 
Therefore, direct photon production in hadron-nucleus scattering provides an excellent test 
of initial state multiple scatterings, while jet or single particle production in photon-nucleus 
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scattering provide independent tests of final state multiple scatterings. Jet and single particle 
production in hadron-nucleus collisions, on the other hand, receive contributions from both 
initial and final state multiple scatterings. Final state multiple scattering in photoproduction 
has been discussed in Ref. 0. Our work will provide the complimentary information on the 
initial state multiple scattering. 

We begin in Sec. |l| with an outline of the formalism used in our calculation. Complete 
analytical results of our calculation are also presented in Sec. ||. The detailed derivation 



of our results, and the calculation of the partonic hard parts are presented in Sec. |T|. In 
Sec. we present our numerical results. We also compare our numerical results with recent 
experimental data. We conclude with a brief summary, and suggestions for further work. 



II. FORMALISM AND ANALYTICAL RESULTS 

An energetic photon can be directly produced at short- distance in high energy collisions, 
and does not interact strongly once produced. Therefore, it has been recognized for a long 
time that direct photon production is a clean probe for short-distance dynamics in high 
energy collisions . Data from hadronic prompt photon production plays a very important 
role in QCD global analysis, and provides constraints on gluon distributions in hadrons 
||8HlO|. In this section, we outline the general formulas for cross sections of direct photon 
production, and present our analytical results for the contribution from double scattering 
processes. 

A. Formalism 

In the following discussion, we study direct photon production in hadron-nucleus colli- 
sions, 

h{p') + A{p) — . 7(0 +X , (1) 
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where p is defined as the averaged momentum per nucleon. In general, the total cross section 
for the above process can be expressed as a sum of contributions from single scattering, 
double scattering and even higher multiple scattering, 

da,A^,{l) = dai'2^^{l) + + ... , (2) 

where the superscripts (S) and (D) represent the single and double scattering, respectively, 
and represents other possible multiple scatterings. In this paper, we consider only the 
double scattering, and its contribution to nuclear dependence. 



As a result of perturbative factorization |Tl|, the single scattering cross section can be 
expressed as 

= AJ2 J dx' fa/hix') J dx fb/wix) daab^yix',xj) . (3) 

a,b 

In Eq. (H), fa/h{x') is a normal parton distributions in the beam hadron h, fb/N{x) is an 
effective nucleon parton distributions inside a nucleus, which should include the well-known 
EMC effect. In principle, the parton-parton scattering cross section, daab^-y, should include 
both direct and fragmentation contributions. That is, an energetic photon can be produced 
directly at short-distances, or produced from the fragmentation of an energetic parton which 
was produced at short-distances []7|,p!2[|. For example, the partonic scattering may produce 
an energetic quark, which radiates a photon. Since we are most interested in fixed target ex- 
periments here, the fragmentation contribution is much smaller than the direct contribution 
in most of phase space |]12[. Therefore, in the rest of our discussion, we will consider only the 
direct production of photons. For example, at the lowest order, we have contributions from 
qq —>■ 7(7 "Annihilation" diagrams, sketched in Fig. |I]a; and gq{oT q) 7g(or q) "Compton" 
diagrams, sketched in Fig. |l|b. 

In terms of the generalized factorization theorem |jl3| , the double scattering cross section 
can be written as: 

d<ji%,{l) = E / dx'fa/H{x') /) , (4) 



where da^^_^^{x' , p, I) can be thought as the double scattering cross section between a parton 
and the nucleus. At the lowest order, it can be factorized as 

da[^^^^{x',p,l) = / dxdxkdxk'J2'^{i}(^^^k,Xk')H{iy{x',x,Xk,Xk',l) • (5) 

{i} 

In Eq. (|^), T^ij{x,Xk,Xk') are the matrix elements of four-parton operators, characterized 
by the set of fields operators {i}; and H^ij are the corresponding partonic hard scattering 
functions. The x, Xk and Xk' are independent coUinear momentum fractions carried by the 
partons from the nucleus. The graphical representation of Eq. (^) is shown in Fig. |^. At the 
lowest order, there are three types of partonic subprocesses that contribute to the double 
scatterings. Feynman diagrams of these partonic subprocesses are sketched in Fig. ^. 

We define the invariant direct photon cross section in hadron-nucleus collision in terms 
of cross sections in hadron-nucleon collisions, 

dH = ^ dH 

dH dH ' 

where Eq. (||) was used. Substituting Eq. into Eq. (j^), we obtain the definition for the 
nuclear dependence parameter 



1 / 1 F ^4aL^(0 

From Eq. (^, a{l) > 1 if dal^_^^/d'^l is positive, which will turn out to be the case for the 
kinematic regime in which we are interested here. However, in general, the double scattering 
contribution a^^^ may be negative, and a{l) < 1 in certain part of phase space. The 
positivity of a cross section requires the sum of all possible multiple scattering contribution 
to be positive. The separation between single and double scattering is not unique. For 
example, two scatterings can be very close to each other, and localized in one nucleon, and 
such a double scattering will not provide the anomalous nuclear dependence, and may be 
classified as a single scattering. 
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We will argue later that the leading double scattering contribution, da^^i^^/d^/, is pro- 
portional to A^^^. Consequently, the value of a{l) will be between 1 and 4/3, depending 
on the relative sizes of contributions from the single and double scatterings. If the double 
scattering contribution is larger than the single scattering contribution in a certain part of 
the phase space, the value of «(/) in that region will approach 4/3. 



B. Analytic Results 



In this subsection, we present the analytic results which are used to calculate the nuclear 
dependence parameter a{l) defined in Eq. (|^. 

Following Eq. the lowest order invariant cross section for single scattering direct 
photon production is given by 



Ej 



da 



hN^-y 



(0 



I dx' fa/h{x') I dxfb/N{x)S ix ^ 

, J J \ X s + u 



dH 



X ttem a. I - 



1 



M. 



(8) 



,s/ \x's + u, 

where X]a,6 I'un over all gluon, quark and antiquark flavors; and the matrix elements for the 
"Annihilation" and "Compton" subprocesses, sketched in Fig. ^ are given by 



2 , /4 



M. 



qg^-yq 



U t \ 



(9a) 
(9b) 



where is the fractional charge carried by a quark of type "g". The invariants s, i and 
u are usual Mandelstam invariants for the parton-parton subprocess. They are related to 
those at the hadron-nucleon interaction by 



s = x'xs, i = x' t, u = xu; 

s = {p' + pf, t={p' ~ If, u = {p~ If 



(10a) 
(10b) 



In the case of double scattering, there are four physical partons linking the matrix ele- 
ments T and the partonic hard parts if, as shown in Eq. (^. After taking into account of 



momentum conservation, there are still three independent momentum fraction integrations 
[x, Xk and Xk' defined in Fig. ^ between the matrix elements and the partonic parts, in 
contrast to one independent momentum fraction integration for the case of single scattering. 
As explained in next section, we take the leading pole approximation to integrate over two 
of the three momentum fractions (e.g., Xk and Xk')- Then, the invariant double scattering 
cross section Eidcr'^^^/d'^l can be reduced into a form very similar to the single scattering 
cross section defined in Eq. (|]). Following our derivation in next section, we obtain 



da 



(D) 



(0 



dH 



aem(47ras)^ / dx'dxS ix 



-x't 



x's + 

E [ fq/h{x') X\ A) Hgq 

+ fq/h{x') A) Hg 

+ fg/h{x')<l>g{X,x',A)Hg 



x's + u 



(11) 



where J2q runs over all quark and antiquark fiavors. In Eq. (|TT]), the functions $4 with 
i = q,g represent the effective parton flux from the nucleus. They are given by 



dx^ 



X 



(x's + u) 



+ 



d mx.AY 
dx \ X 



—u 



x' six's + u) 



:i2) 



The Ti{x,A) with i = q, g in Eq. (|I^) are the twist-four matrix elements in nuclei. They 
were originally introduced in Ref. [^, and are given by 



T,(x,A) 



2tx 



27r 



(13a) 



and 



Tg{X,A) 



X ■ 



<kh_ 
1 

XJ)^ 



dy__dy2 

2n 



■0{yi-y2n-y 



{pA\F^^{y:,)F^^{Q)F^^{y^)FX{y-)\pA) 



(13b) 



In Eq. ([131), and ipq are the fleld strength and quark fleld operator, respectively. 
In Eq. (p!TD, the Hi are the partonic hard parts, and 



/ 2 \ / -u x's + u\ . ^ , 

^^,= (41 3^ + ^) ; (14b) 



,16/ \ — M x's 
which we will derive in next section 



36/ \x's + u 

+ (14c) 



Eqs. (|TT]), (|12D and ([l^) are our complete analytic results at leading nonvanishing order 
in Ug- As usual, the next-to-leading order (NLO) contribution might be important for the 
single and/or double scatterings. However, since the nuclear dependence parameter a{l), 
defined in Eq. (|^), depends on the ratio of the double and single scattering contributions, 
we expect that the values of a{l) presented in this paper are not very sensitive to the NLO 
contributions. 



III. DERIVATION OF THE DOUBLE SCATTERING CONTRIBUTIONS 

In this section we provide the derivation that leads to the analytic results presented in last 
section. The method that we used here was first introduced in Ref. 0. It can be summarized 
in following technical steps: a) factorize the double scattering contribution into a convolution 
between the partonic hard parts and the corresponding multiparton matrix elements (e.g., 
see Eq. (H)); b) in the leading pole approximation, integrate over two of the three independent 
momentum fractions by contour integrations, and reexpress the multiparton matrix elements 
in terms of the Tq{x,A) and Tg{x,A) defined in Eq. (p!3|); c) calculate the corresponding 
partonic hard parts. 

At lowest order, only three types of partonic subprocesses, as sketched in Fig. ^ con- 
tribute to the double scattering cross section da'f^^^ introduced in Eq. (^. These three sub- 
processes correspond to adding two gluons to the lowest order "Annihilation" and "Comp- 
ton" subprocesses, shown in Fig. |T]. In the following subsections, we present the detailed 
derivation for one subprocess, and provide the results for other subprocesses. 
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A. Perturbative factorization 



Consider the subprocess shown in Fig. in which there are three independent 
four- momentum hnking the partonic part and corresponding two-quark-two-gluon matrix 
element. In the center of mass frame of high energy colhsion, all partons inside the nucleus 
are moving almost parallel to each other, along the direction of the nucleus. Therefore, all 
three parton momenta can be approximately replaced by the components coUinear to the 
hadron momentum. After such a coUinear expansion, the double scattering contribution 
from the generalized "Annihilation" subprocess shown in Fig. |^a can be written as 

El ^3p^ = j dxdxkdxk' J d'^kTT{x,Xk,Xk',kT,p) H{x'p',x,Xk,Xk',kT,p,l) , (15) 

where 2x's is the flux factor between the incoming beam quark and the nucleus, and x'p' 
is the momentum carried by the beam quark. In Eq. (p!5D , the two-quark-two-gluon matrix 
element, T, is defined as 



T{x,Xk,Xk',kT,p) 



dVi dy2 dy d^yr 
27r 27r 27r (27r)2 

x^{PA\A+iy2,OT)Mo)i^Myi)A-^iy'^yT)\PA) ■ (i6) 

The corresponding partonic part H is given by the diagrams shown in Fig. ^, with gluon 
lines contracted with ppp'^, quark lines from the target traced with {'y ■p)/2, and quark lines 
from the beam traced with (7 ■ {x'p'))/2. Here, we work in Feynman gauge, in which the 
leading contribution from the gluon field operators is A'' ^ A'^ {p'' / p'^) . We also kept the /ct 
for the gluons in order to extract a double scattering contribution beyond leading twist. 
By expanding the partonic part H introduced in Eq. ([I5| ) at /ct = 0, we have 



H{x'p',x,Xk,Xk',kT,pJ) = H{x'p',x,Xk,Xk',kT = 0,p, /) 

1 dm 



dH 

+ 



dk^n 



krp -\- 



kT=o ^ dkj.dk^ 



k^ki^ + ... . (17) 

kT=0 
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In the right-hand-side of Eq. ([TTD, the first term is the leading twist eikonal contribution, 
which does not correspond to physical double scattering, but simply makes the single- 
scattering matrix element gauge invariant. The second term vanishes after integrating over 
/ct- The third term will give the finite contribution to the multiple scattering process. 
Substituting Eq. ([T^ into Eq. and integrating over (Phx, we obtain 



E, 



da 



dH 



2x's 
1 



X 



-2' 



a/3 



dx dxk dxk' T{x, Xki x^', A) 
1 



2 dk^dk^ 



H{x'p',x,Xk,Xk', kr = 0,p, /) 



(18) 



where the modified matrix element T is given by 

T{x,Xk,Xk',A) = 



dyi_ dy_ dy2_ ^i^p+y- ^ix^v+y- ^-i{xf,-x,^,)p+y- 

27r 271 271 



x-{p^\F^+{y,)ij,{0)^+^,{y-)F+-{y-)\pA) 



(19) 



In Eq. (p^), = F^'^nf^, and F^" is the field strength, and vector rijs = 6(^+. 

In following sections, we show how to perform the integrations over parton momentum 
fractions, and evaluate the partonic parts {d^ / dk^dkj,) H for different subprocesses. 



B. Leading pole approximation 

The double scattering contribution defined in Eq. (|18D depends on integrations over 
three partonic momentum fractions x, x^, Xy- If all partons in Fig. ^ carry some finite 
momentum fractions, the oscillations of the exponentials in the matrix element T defined in 
Eq. ([1^) will destroy any nuclear size enhancement that could come from the y integrations. 
However, even at the lowest order, we find that there are some Feynman diagrams which 
have two poles corresponding to zero momentum fraction partons, and, these poles are not 
pinched. Therefore, two of the three parton momentum fractions can be integrated explicitly 
by contour integration. These integrations will eliminate two exponentials, and thus, the 
corresponding y integration could provide the nuclear size enhancement up to A^^^ . But, in 
terms of double scattering picture, if we require two soft field operators to come from the 
same nucleon, we will get the familiar A^^^ enhancement. 
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Of course, there are double scattering diagrams without such poles, but we expect A" 
(a > 1) dependence only when the poles are present. In this paper, we evaluate only 
diagrams that have such poles, and we call our results at leading pole approximation. 

In order to perform the integration of momentum fractions, it is convenient to rewrite 
the double scattering contribution defined in Eq. (|T^) as 

dalj^l^ 1 r dyi dy' dy2 1 



(20) 



X I -is"" 



In Eq. (|20|) , the modified partonic part H is defined as 



H{yi,y-,y2,kT,P,l) = J dxdXkdxy e'^^^^y^ ^Ix^p+y- ^-i{x,-xy)p+y^ 

xH{x'p',x,Xk,Xk',kT,p,l) , (21) 



where the partonic part H is given by diagrams shown in Fig. ^ It is clear from Eq. (pT]) 
that all integrals of momentum fractions can now be done explicitly without knowing the 
details of the multiparton matrix elements. 

Consider the diagram shown in Fig. The final state photon-gluon two particle phase 
space can be written as 

^11/ x't -k^-2kT-l\ 

r = — 6\x + Xk + + — — . (22) 

Svr"^ x's + u \ x's + u x's + u J 

In deriving Eq. (^), we have omitted the factor d^l/Ei, due to the definition of the invariant 
cross section (e.g., see Eq. ([20|) ). Using Eq. (|2^), the contribution to H from the diagram 
shown in Fig. ^ can be expressed as 

Hi -a = — ■ Hi_a{x, Xk, Xk') 

2tt x's + u 



1 1 x't ~B-2kT-V 
fc^ P d X + Xfc + ^ + — 

Xk - Xk' - J; - Xk - + le \ xs + u x s + u ^ 



where the subscript "/-a" has following convention: "/" stands for the type-/ subprocess, 
shown in Fig. "a" for the real contribution, corresponding to diagrams in Fig. In 
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Eq. (p3D, the factor C/ is an overall color factor for the type-/ subprocess. The function 
Hi-a in Eq. ( pB]) is given by 
1 1 



Hi-a = - — Tr 
4 x's 



7 ■ {x'p + A;t)7 ■ 1^7 ■ {^'p' + kT)Ri'^a7 -pLj^a {-9ap) {-9iiu) , (24) 



where -R/^^q and /v/^a the right and left blob, respectively, as shown in Fig. |a. These 
blobs include all possible tree Feynman diagrams with the external partons shown in the 
figure. Substituting Eq. ( p3D into Eq. (^), we obtain 

1 r , 1 /• . 1 



2tt x's + u 



Xk-§; + te 



Xk - Xk' - ^ - ie 



X / dxe'^'P y^ 5 [x + Xk + 



x't 



+ 



-kl-2kT-V 



Hl-a{x,Xk,Xk' 



(25) 



x's + u x's + u 

After performing dxk and (ia;^' by contour integration, and dx by the 5-function, we derive 



1 



x6'(-y2") 6'(?/j; - y") Hl-aix, Xk, Xk') , 



i//_a = (27ra,) , , e-- e 



(26) 



where the 6'-functions result from the contour integrations, and the momentum fractions for 
the function Hi-a are defined as 

1 



X 



x's + u 



u 

x's 



xt+ -^k^ -2kT-l 



kj^ 



Xk , ) 

x's 



Xk' = ; 



x't 



X 



x's + u 

Similarly, we derive contribution from the diagram shown in Fig. ^ as 

1 



(27a) 

(27b) 
(27c) 
(27d) 



i^V^Vl g«(fcTA's)p+(y -J/2 ) 

^G{.y2 - y') ^ivi - V2) Hi-b{,x, Xk, Xk') , 



Hi_b = (27ra,) Ci -— — e'^'' e 
x's + u 



(28) 



where x, Xk and Xk' are also defined in Eq. (P?]). Similarly to Eq. (P^, the partonic part 
Hi^b is given by 
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7 ■ (xV)i?^^b7 ■ {-9ap) {-9t.u) ■ 



(29) 



The diagram shown in Fig. ^ has following contribution 



1 



xe{y~ -y^)9{-y')Hj_b{x,Xk,Xk') • 



Hi_^ = (27ra,) e^'' e 

x's + u 



(30) 



In deriving Eq. (0), we used the fact that the partonic part Hi-c = Hj_i, when Xk and x^/ 
are evaluated at the same values as listed in Eq. (pl). 



Combining Hj^a, Hj^b and Hj^c (given in Eqs. (|26|) , (^), and (|30D , respectively) to- 
gether, we obtain the total contribution to H, defined in Eq. (|2T|), from the type-/ diagrams 
shown in Fig. 

1 



X 



(31) 



All momentum fractions in Eq. (|3l|) are evaluated at the values defined in Eq. (pTj). In 
deriving Eq. (|3lD , we have dropped a term proportional to 



[0{-y2) 0{y, - y-) - e{y, - y-) d{y^ - y,) - 9{y- - y,) 9{-y-) ] 



. 



This is because of the phase exp[ixp'^y~] which effectively restricts ~ l/{xp^) 0. 
Physically, it means that all y integrations in such term are localized, and therefore, will not 
give any large nuclear size enhancement. 

By substituting Eq. (|3T|) into Eq. (|20|), we can obtain the lowest order double scattering 
contribution from the type-/ diagrams shown in Fig. One important step in getting the 
final result is taking the derivative with respect to /ct as defined in Eq. (^). Comparing 
Eq. ( pTD with Eq. (0), and observing that 



^iSp+y, Hi_aix,Xk,Xk') - e'^'P^y^ Hl-aix,Xk,Xk' 



kT=0 







(32) 
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we found that the derivatives on the exponential exp[i {k'^ / x' s) p'^ {y~ — y^)] do not con- 
tribute, and that we can therefore set exp[z {k^ / x' s) {y^ — ^2^)] = 1 in Eq. Substi- 
tuting Eq. ( |3TD into Eq. (^), and use Eq. ( |13a| ), we obtain 

El = aem[^nas) e 6/ 



(PI 2x's x's + u 

1 «\ 1 ^2 _ .. _ 

Tq{x, A) Hi^a{x, Xk, Xk') - Tq{x, A) Hi_h{x, Xk, Xk') , (33) 



X I — ^"'^ 



2 J 2 dk^dk^ 

where a\^^ stands for the double scattering contribution from the type-/ subprocess shown 
in Fig. It is important to note that although the interference diagrams shown in Fig. ^ 
and Fig. ^c are important in driving Eq. (^), the final result depends only on the real 
diagram shown in Fig. That is, the double scattering picture is preserved. The role of 
interference diagrams is to take care of the infrared sensitivities of the short-distance hard 
parts. 



C. Final factorized form 



The derivatives with respect to kr in Eq. (|33|) are straightforward. It is most convenient 
to reexpress derivatives with respect to kx in terms of derivatives with respect to x or x. 
After working out the derivatives, we obtain H 



da 



fTq{x,Ay 



x2 



dx'^ 



(x's + u) 



+ 



d_ mx,A) - 

dx \ X 



—u 



x' six's + u) 



where x is given in Eq. (|27d| ), and where the partonic hard part Hqq is defined as 

Hqq = ClX Hl-a{x, Xk = 0, Xk' = 0) . 



(34) 



(35) 



Following the same derivation, we obtain contributions from the type-// and tjpe-III 
diagrams shown in Fig. ^. For the type-// diagrams, as sketched in Fig. Qd, we have 



da 



(D) 
II 



aem(47ra, 



\2 „2 



x2 



d^l 2x's x's + u 

fTq{x,Ay 



dx"^ 



X 



(x's + uY 



+ 



d ( T,{x,Ar 
dx \ 



X 



-u 



x' six's + u) 



(36) 



15 



In Eq. (|36|), the partonic hard part Hg is defined as 



Hg = Cii Hn-a{x, Xk = 0, Xk' = 0) , 



(37) 



where Cu is the overall color factor for the type-// diagrams, and Hn-a is given by the real 
diagrams shown in Fig. ^, and defined as 



1 



Hii-a = 7 Tr 7 ■ {x'p' + XkP + kT)Ru-al " ^' ^//-a {-9ap) {-9t,u) , 



(38) 



where /' = x'p' + (x + Xk)p — I is the momentum carried by the quark going to final state. 
Similarly, for the tjpe-III diagrams, as sketched in Fig. we obtain 



El 



da 



III 



x2 



«em(4vras) ^^^^ 
(T,ix,A) 



^//, 

2x's x's + M 



dx'^ 



X 



_ (x's + uY 



d_ f T,{x,Ay 
dx \ 



X 



-u 



x's{x's + u) ^ 



(39) 



The partonic hard part Hg in Eq. (^) is defined as 



Hq = CinxHni-a{x,Xk = 0,Xk' = 0) 



(40) 



where Cm is the overall color factor for the type-3 diagrams, and Hm-a is given by the 
real diagrams shown in Fig. ^, and defined as 

1 



Hiii-a = ^ Tr 



7 -pRui-al ■ I' L'^III-a i-9af3) (-^^^) 



(41) 



where /' is the same as that defined after Eq. (|38D . 

The partonic short-distance hard parts, defined in Eqs. (|35|) , (|37D and (|40|), can be easily 
evaluated by calculating the corresponding Feynman diagrams. Our results were presented 



in Eq. (|T|) in Sec. |ITB| . 

After convoluting Eqs. (p4D, (|36D, and (|39D with the corresponding parton distributions 
from the beam, we obtain the complete analytical expressions for the double scattering 
contribution in hadron-nucleus collisions, which was presented in Eq. (|Tl]) of Sec. [11 B . 
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IV. NUMERICAL RESULTS AND DISCUSSIONS 



In this section, we present our numerical results for the Cronin effect in direct photon 
production. We numerically evaluate the nuclear dependence parameter a{l) defined in 
Eq. by using our analytical results presented in Eqs. (H) and (0)) we also compare 
our numerical results with recent data from Fermilab experiment E706 . 

The nuclear dependence parameter a{l) defined in Eq. depends on contributions from 
both single scattering and double scattering. All these contributions depend on the nonper- 
turbative parton distributions or multi-parton correlation functions. In deriving following 



numerical results, the Set 1 pion distributions of Ref. |TJ] are used for pion beams; and the 
CTEQ3L parton distributions of Ref. are used for free nucleons. The twist-4 multi-parton 
correlation functions defined in Eq. ( [T3| ) have not been well- measured yet. By comparing 
the definition of these twist-4 correlation functions with the normal twist-2 parton distri- 
butions |JT5|], authors of Ref. ||^ proposed following approximate expressions for the twist-4 



correlation functions, 

Ti{x,A) = \^A'/'f,/Aix,A) (42) 

where i = q,q, and g. The fi/A are the effective twist-2 parton distributions in nuclei, and 
the factor A^^^ is proportional to the size of nucleus. The constant has dimensions of 
[energy]^ due to the difference between twist-4 and twist-2 matrix elements. The value of 
was estimated in Ref. |[16| by using the measured nuclear enhancement of the momentum 



imbalance of two jets in photon-nucleus collisions [|1^,|T^, and was found to be 



A^ ~ 0.05 - O.lGeV^ . (43) 

This value is not too far away from the naive expectation from the dimensional analysis, 
A^ ~ ^QCD- calculation below, we use A^ = 0.1 GeV^. Therefore, our numerical 

results can be thought as the upper limit of the theoretical predictions. 

The A^^^ dependence of the twist-4 multiparton correlation functions, introduced in 
Eq. (|4^), is not unique. From the definition of the correlation functions in Eq. ([T3|) , the lack 
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of oscillation factors for both y~ and integrals can in principle give nuclear enhancement 
proportional to A^/'^. The A^^^ dependence is a result of the assumption that the positions 
of two field strengths (at and respectively) are confined within one nucleon. 

In Eq. (^2]), the effective nuclear parton distributions fi/A should have the same operator 
definitions of the normal parton distributions with free nucleon states replaced by the nuclear 
states. For a nucleus with Z protons and atomic number A, we define 

U/a{x, A) = A {^U/n{x) + |/./p(x)) i?™^(a;, A) , (44) 

where fi/N{x) and fi/p{x) with i = q,q,g are normal parton distributions in a free neutron 
and proton, respectively; and N = A — Z . The factor Rf^'^ takes care of the EMC effect in 
these effective nuclear parton distributions. We adopted the fj^Qj^i Ref. [|19|, which fits 



the data well. However, at fixed target energies, the x values covered by the direct photon 
experiments are large and out of the nuclear shadowing region. The integration over dx' in 
Eqs. (^ and (^) averages out the EMC effect from the large x region. Actually, one can 
neglect the Rf^^ in Eq. (||). 

Using the parton distributions and correlation functions introduced above, and our ana- 
lytic results presented in Eqs. (|^) and (|TT]), we can derive the nuclear dependence parameter 
a{l), defined in Eq. (|^), without any further free parameter. 

In Fig. ^ we compare our numerical predictions for the nuclear dependence parame- 
ter with the recent experimental data from Fermilab experiment E706 0. The aE706(0 
presented in Fig. |^ is slightly different from that defined in Eq. (^). E706 measured the 
direct photon cross sections with the vr" beam on two different targets: Cu(A = 63.55) and 
Be{A = 9.01); and the aE706(0 extracted according to following definition 

The beam energy is p' = 515 GeV. It is clear that the theoretical calculation presented in 
this work is consistent with the data. 

It is evident from Fig. |^ that the nuclear dependence parameter aE706(0 is very close to 
the unity, or equivalently, the Cronin effect for direct photon production is very small, and 



much smaller than that observed in the single particle inclusive cross sections |T[ . One clear 
difference between the direct photon and the single particle inclusive cross sections is that 
direct photon production has only initial state multiple scattering, while the single particle 
inclusive has both initial and final state multiple scattering. In addition, the single particle 
inclusive cross sections depend on the parton-to-hadron fragmentation functions. 

As pointed out in Ref. P], the multiple scattering contribution is most important when 
the momentum fraction x from the nuclear correlation functions is large because of the 
derivatives with respective to the x, which were introduced in Eq. ([I^). However, for direct 
photon production, the kinematics do not fix all parton momentum fractions, and leave one 
momentum fraction to be integrated, for example, x' in Eqs. (|]) and (pHI). Because of the 
steeply falling nature of the distributions and correlation functions, the cross sections in 
the central region {xp ^ 0) for a given value of It is dominated by the distributions with 
momentum fractions x' ~ x ~ xt = 2/^/-^/?, which is less than 0.6 even for the largest 
value of It shown in Fig. |^. Therefore, the double scattering contribution is relatively small 
because the derivative terms are not significantly enhanced, and consequently, aE706(0 
close to one. 

In contrast, for inclusive single hadron production, the parton-to-hadron fragmentation 
functions effectively shift the contribution at a given It to the large x region because of all 
fragmentation functions vanish when z goes to 1. Kinematically, direct photon production 
corresponds to single particle production at 2; = 1. We, therefore, expect that single hadron 
production has a larger Cronin effect than the direct photon production at the exact same 
kinematics, even before including contributions of final state multiple scattering. 

In the case of a pion beam, the quark-antiquark "annihilation" subprocess, sketched in 
Fig. |l|a, dominates the production of direct photons at the fixed target energy, due to the 
valence antiquarks in the beam. However, for a proton beam, the quark-gluon "Compton" 
subprocess, as sketched in Fig. ^d, is more important for the production of direct photons. 
Therefore, direct photon production with a proton beam is more sensitive to the proton 
gluon distributions. In Figs. ||, and we present our predictions of nuclear dependence 
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parameter a{l), defined in Eq. (|^), for a vr" and proton beam, respectively. 

In Figs. ||a and ||a, tlie nuclear dependence parameter a{l), defined in Eq. (|^, is plotted 
as a function of photon's transverse momentum It at xp = 0. In plotting these figures, a 515 
GeV beam energy and a Copper target was assumed, using the same parton distributions 
and correlation functions as in Fig. |^. As expected, and as found in Fig. |^, the value of 
a{l) is very close to one for both pion and proton beams. Changing a pion beam by a 
proton beam does not affect the kinematics of the collisions. The effective values of parton 
momentum fractions from the nuclear target are the same for both cases. Therefore, as 
explained above, the values of a{l) is close to one due to the fact that the effective parton 
momentum fractions from the nuclear target are not large. 

In order to enhance the contribution from the double scattering, we need to look for 
events at large negative xp, where the effective values of parton momentum fractions from 
the target are larger. In Figs. ||b and ^3, we plot the nuclear dependence parameter a{l), 
defined in Eq. (^, as a function of xp at It = 6.0 GeV. The same beam, target and beam 
energy were used. It is clear that when xp becomes large and negative, a{l) increases. This 
is consistent with the fact that the negatively larger values of xp, the bigger effective values 
of parton momentum fractions from the nuclear target, and consequently, larger derivatives, 
defined in Eq. (|12D. 

Comparing Figs. Pd and Pd, one finds that as xp decreases, the values of a{l) with a 
proton beam increases much fast than that with a pion beam. This is because the Compton 
subprocess dominates the production of direct photons in the case of a proton beam, and the 
gluon distribution in a proton falls much faster than the valence quark distributions as the 
momentum fraction increases. The more rapidly falling gluon distribution produces larger 
derivative terms, and therefore, larger values of a{l). Future data from Fermilab experiment 
E706 with a proton beam can test this feature. 

In summary, using generalized factorization in QCD perturbation theory, and using the 
method, developed in Ref. for calculating nuclear enhancements, we demonstrated in this 
paper that the observed small Cronin effect in direct photon production is consistent with the 
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much larger Cronin effect observed in single jet and single particle inclusive cross sections. 
We hope that the same method can be used to explain the puzzle for the nuclear dependence 
in the momentum imbalance. Data on dijet momentum imbalance in photoproduction 
and hadroproduction [|I^ have shown strong nuclear dependence, while much smaller nuclear 
dependence has been seen in momentum imbalance of Drell-Yan pairs PD|. 
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FIGURES 

FIG. 1. Lowest order Feynman diagrams contribute to single scattering: a) "Annihilation", b) 
"Compton". 

FIG. 2. A Graphical representation of double scattering contributions from the parton-nucleus 
collisions. 

FIG. 3. Three types of leading order Feynman diagrams contribute to the double scattering. 

a) Type-I: "Annihilation" diagrams corresponding to the two-quark-two-gluon matrix element; 

b) Type-II: "Compton" diagrams corresponding to the four-gluon matrix element; c) Type-Ill: 
"Compton" diagrams corresponding to the two-quark-two-gluon matrix element. 

FIG. 4. Feynman diagrams for the "Annihilation" diagrams corresponding to the 
two-quark-two-gluon matrix element: a) real diagrams, b) and c) interference diagrams. 

FIG. 5. The real "Compton" diagrams corresponding to the four-gluon matrix element. 

FIG. 6. The real "Compton" diagrams corresponding to the two-quark-two-gluon matrix ele- 
ment. 

FIG. 7. Cronin effect in direct photon production with 515 GeV vr" beam on Cu and Be targets. 
The theory curve is from Eq. (^), and the data is from Fermilab experiment E706 

FIG. 8. Theoretical predictions of a{l) of Eq. function of photon's transverse 

momentum at xp = 0.0; b) as a function oi xp at It = 6.0 GeV. 515 GeV vr^ beam and Copper 
target were used. 

FIG. 9. Theoretical predictions of a{l) of Eq. (^), as in Fig. |8|, but with a proton beam. 
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